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Time-dependent wave-packet approach for fusion reactions of halo nuclei 
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The fusion reaction of a halo nucleus n Be on 208 Pb is described by a three-body direct 
reaction model. A time-dependent wave packet approach is applied to a three-body 
reaction problem. The wave packet approach enables us to obtain scattering solutions 
without considering the three-body scattering boundary conditions. The time evolution 
of the wave packet also helps us to obtain intuitive understanding of the reaction dynamics. 
q ; The calculations indicate a decrease of the fusion probability by the presence of the halo 
CO ■ neutron. 
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1. INTRODUCTION 

Since the discovery of the halo nuclei, much effort has been devoted to develop reaction 
theories appropriate for weakly bound projectiles. The eikonal approximation has been 
found to be useful for reactions at medium and high incident energies lj. However, even 
a basic picture has not yet been established for the reactions at low incident energies. 
For example, there have been controversial arguments on whether the fusion probability 
is enhanced or suppressed by the presence of the halo nucleon. 

We have been analyzing the low energy reactions of halo nuclei in a three-body direct 
reaction model. A time-dependent wave packet approach has been developed for this 
purpose. We first analyzed a one-dimensional three-body model [2]. We then made anal- 
ysis of the fusion reaction of n Be on a medium mass target, 40 Ca, in a three-dimensional 
model [3 . The calculation was restricted to the cases of total angular momentum being 
equal to zero (head-on collision). In these analyses, we have found that the fusion proba- 
bility is hindered by the presence of the halo neutron. This conclusion is opposite to many 
other theoretical works which have claimed the increase of the fusion probability |4l5j . 

We here report our extended analyses with the wave packet method for the reaction of 
n Be on 208 Pb. For this reaction, the analyses have been made recently with the three- 
body model by other groups [H|T] . The measurements are also available for systems close 
to thisjE]. 

The Hamiltonian of the three-body model is time-independent. Still, there are several 
reasons to employ the time-dependent approach for the static problem. The time evolution 
of the wave packet gives us intuitive pictures for the reaction mechanism. It is not 
necessary to prepare complicated scattering boundary conditions for the three-body final 
states. The reaction probabilities for a certain energy region can be obtained at once from 
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a single wave-packet solution. The trade-off for these advantages is a heavy computational 
cost. The method has also been developed and applied in the field of chemical reactions[9j. 

2. THE THREE-BODY MODEL 

We describe a reaction of a single-neutron halo nucleus by the three-body model con- 
sisting of the halo neutron (n), the core nucleus (C), and the target nucleus (T). The 
projectile (P) is composed of the halo neutron and the core nucleus, P=n+C. We con- 
sider the reaction of n Be on 208 Pb. Namely, the core is the 10 Be and the target is the 
208 Pb. The time-dependent Schrodinger equation is expressed as 
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where we denote the relative n-C coordinate as r and the relative P-T coordinate as 
R. The reduced masses of n-C and P-T motions are m and fi, respectively. The n-C 
potential, V n c(r), is real. The n-T potential, V n T(r nT ), is also chosen to be real. The C-T 
potential, Vct{Rct), is complex. The real part of Vct includes the Coulomb and nuclear 
potentials. The imaginary part of Vct describes the fusion process as a loss of the flux in 
the reaction between the core and the target nuclei. All the nuclear potentials are taken 
to be of Woods-Saxon shape. For the n-C potential, r nC = 1.3 fm and a nC = 0.7 fm. 
The strength of the potential is set to give the 2s orbital energy at —0.6 MeV, close to 
the neutron separation energy of the n Be nucleus. For the n-T potential, r n T = 1.27 fm 
and a n T = 0.67 fm. The strength will be varied. For the real part of the C-T potential, 
we employ Tq t = 1.192 fm, ag T = 0.63 fm, and V CT = -46.764 MeV. The Coulomb 
potential is that of the uniformly charged sphere with the radius parameter, r^ T = 1.26 
fm. For the imaginary part of the C-T potential, we employ r CT = 0.8 fm, a CT = 0.4 fm, 
and W C t = -50 MeV. 

In this three-body model, we define the fusion process of the three-body reaction as 
a loss of the flux caused by the absorptive C-T potential. Namely, the fusion is sup- 
posed to occur between the core and the target nuclei when they come closer beyond the 
Coulomb barrier. In this definition, the fusion probability includes both the complete and 
incomplete fusions, since the lost flux does not distinguish final states of the neutron. 

In practical calculations, we make a partial wave expansion. There are two options 
for this; the body-fixed representation and the space-fixed one|10j. For computations at 
finite total angular momenta, the former is superior to the latter. Here we only consider 
the reactions with zero total angular momentum, for which both representations merge 
into the same expression, 

ij) (R,r,t)=>^ — ——^(cosfl), (2) 

where 9 represents the angle between the two vectors R and r. The I is the relative 
angular momentum of the n-C motion, which is now identical to that of the P-T motion 
in the case of J = 0. In numerical calculations, one has to restrict the summation within 
a finite range, I < l ma x- 
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As the initial wave packet, we employ a Gaussian packet for the P-T motion, 

Uq(R, r, t ) = exp[- 7 ( J R - R ) 2 - iK Q R](j) {r) , (3) 

where 0o( r ) is the 2s orbital of the n-C motion and ui = for I ^ at t = t . The Rq 
specifies the P-T separation at the center of the wave packet and K Q indicates the P-T 
wave number on average at t — t^. This wave packet state has the average total energy 
approximately given by E = Ti 2 Kl/2p + Vct(Ro) — e o? where e is the binding energy of 
the orbital 4>o- The Rq is set to be large so that Vct and V u t vanish except for the C-T 
Coulomb potential. We set Rq = 25 fm and 7 = 0.05 fm -2 . 

To solve the time-dependent Schrodinger equation for the radial wave function ui(R, r, t), 
we discretize the radial variables. We treat radial region < R < 50 fm and < r < 60 
fm with radial step AR = 0.5 fm and Ar = 1 fm, respectively. The discrete variables 
representation is used for the second-order differential operator [TT]. It is important to 
employ very high partial waves / to describe the transfer reactions in the Jacobi coordi- 
nate of the incident channel. We take up to l max = 70. The wave function m(R,r,t) is 
thus described on the grid points of about 420,000. For the time evolution, we employ 
the Taylor expansion method. Previously we employed the Crank-Nicholson formula[3 . 
The Taylor expansion method requires usage of the smaller time step and gives more ac- 
curate results. The time evolution is continued until most parts of the wave packet leave 
the interaction region after scatterings and the potentials from the target again becomes 
negligible except for the Coulomb. 



3. WAVE PACKET DYNAMICS 

The wave packet is a superposition of the solutions for various total energies around Eq. 
It provides us an intuitive picture for the reaction dynamics around this energy. In Fig. ^ 
we show a time evolution of the wave packet. The K value is so chosen that the average 
total energy E is about 38 MeV, which is close to the barrier top energy. The left-panels 
show the density distribution integrated over the angle 9, p(R, r, t) = J2i \ u i(R, r > t)\ 2 - 

Here the abscissa is r and the ordinate is R. The right-panels show the density distribution 

r 1 V 2 

integrated over R, p(r,9,t) = ^JdR\ ^iUi(R,r,t) ^ P;(cos#)| 2 . Here the abscissa 
is x = r cos 9 and the ordinate is y = r sin 9. The center corresponds to the center-of-mass 
of the projectile, and the x direction is parallel to the R vector. In the figure, the target 
nucleus approaches to the center from the right, then returns back. The time evolves from 
the top to the bottom. 

Let us first look at the five panels in the left. The top- left panel shows the initial wave 
packet. There is a nodal structure of the 2s orbital in the r direction. The distribution 
extends to large r reflecting the halo structure. In the R direction, the wave packet shows 
the Gaussian distribution around Rq. As the time evolves, the wave packet moves towards 
smaller R. The third (middle) panel roughly corresponds to the classical turning point. 
Some components of the wave packet pass through the C-T Coulomb barrier, while some 
are reflected by the barrier. The components going beyond the barrier are absorbed by 
the imaginary C-T potential and disappear. After the wave packet reaches the closest 
point, a flow to larger r is seen. This corresponds to the breakup process. Performing 
the calculation with the n-T potential being switched off, we have confirmed that this 
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Figure 1. The time evolution of the wave packet solution. In the left panels, the density 
is plotted as a function of R and r. The abscissa is the P-T distance R, and the ordinate 
is the n-C distance r. In the right panels, the density is plotted as a function of r and 9 
of n-C coordinate r, after integrating over R. The initial wave packet is shown in the top 
panel, and the final state in the bottom. See text for more detail. 
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is induced by the C-T Coulomb force. One can also see the wave packet component in 
the diagonal (R ~ r) direction. This is a transfer component. By taking account of high 
partial waves (l max = 70), one can describe the transfer processes in the Jacobi coordinate 
system of the incident channel. A nodal structure near the right end in the bottom-left 
panel is due to a spurious reflection at the boundary in the r coordinate. 

The five panels in the right-hand side are useful to see behaviors of the halo neutron. 
The initial distribution is spherical since we assume 2s orbital in n Be. Because of the 
acceleration of the 10 Be core induced by the target Coulomb field, the halo neutron appears 
to be accelerated to the right (towards the target). Note that the center in these figures 
corresponds to the center of mass of the projectile. This acceleration induces the Coulomb 
breakup. As the target nucleus gets closer, one can also see the formation of molecule-like 
orbital in the third (middle) panel. In the bottom panel, the breakup component shows 
the fan-shaped, multiple waves. This is again caused by the spurious reflections at the 
boundary in r coordinate. The transfer component is also seen as an oval distribution 
which locates its center at the average position of the target nucleus. 

4. FUSION PROBABILITY 

To obtain the fusion probability for a fixed incident energy from the wave packet so- 
lution, we apply the energy projection procedure. For a wave packet state ip a , we define 
the energy distribution function W a (E) as 

W a (E) = (ip a \5(E — H)\ip a ), (4) 

where H is the three-body Hamiltonian in Eq. (Q). We evaluate Eq. (0J) for the initial 
(a = i) and the final wave packet (a = f). The fusion probability is then obtained as a 
fraction of the energy distribution that was absorbed by the C-T absorptive potential. 

_ Wj(E) — Wf(E) 

Pfusion(£) - W ^ F) • (5) 

In calculating the energy distribution function W a (E), one may also employ the time 
evolution of the wave packet. Fourier transforming the delta function, one obtains 

^) = ^Ref (~)l*Q)>, (6) 

where ip a if) is the solution of the time-dependent Schrodinger equation (JTJ) with the initial 
condition ^ a (0) = ip a . 

Let us elucidate a basic picture of the reaction dynamics from the calculated fusion 
probabilities. Figure El shows the fusion probability of central collision (zero total angular 
momentum) as a function of the P-T relative incident energy. The fusion probabilities 
of three-body calculations are compared with that of the two-body 10 Be- 208 Pb reaction 
without a halo neutron. The three-body calculations have been carried out with different 
n-T potentials: (i) V n T = (dotted curve), (ii) V n T = —48 MeV (dot-dashed), and (iii) 
V n T = —42 MeV (dot-dot-dashed). The Fig. [T] corresponds to the case (iii), where the 
substantial transfer component is seen. There is small transfer component in the case of 
(ii). 
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Figure 2. Fusion probability SIS db func- 
tion of the projectile-target relative inci- 
dent energy for the case of zero total angu- 
lar momentum. Results of various three- 
body calculations are compared. See text 
for the details. 



Figure 3. Fusion probability func- 
tion of the projectle-target relative inci- 
dent energy. Dependence on the partial 
wave truncation is shown. See text for the 
details. 



As is seen from the figure, the fusion probabilities in the three-body calculation are 
always smaller than that of the two-body calculation. In other words, the fusion proba- 
bility is suppressed by the presence of the halo neutron. Furthermore, this suppression is 
caused principally by the C-T Coulomb interaction. The n-T potential plays a minor role 
in this fusion reaction. 

We could imagine two possible mechanisms that are responsible for the fusion sup- 
pression. The first one is a spectator role of the halo neutron. As is seen in Fig. ^ the 
halo neutron is emitted to the forwared direction. This indicates that the halo neutron 
proceeds straight keeping the incident velocity, while the core nucleus is reflected by the 
target Coulomb field. In this spectator picture, the relative energy between the core and 
the target nuclei is effectively smaller than the projectile-target relative energy. This 
mechanism, which we discussed previously [3], is expected to shift the energy dependence 
of the fusion probability by a constant amount. The calculated fusion probability, how- 
ever, does not show such a simple shift in energy but shows the stronger suppression above 
the barrier. The second mechanism that could explain this property is the excitation of 
the n-C motion. The Coulomb interaction between the core and target nuclei induces 
the internal excitation of the projectile. The excitation energy will distribute in a certain 
energy region with a fluctuation. The energy conservation requires that the P-T relative 
energy decreases to cancel the projectile excitation energy. This results in the decrease 
of the effective incident energy with the fluctuation. The suppression of the fusion prob- 
ability persists even high above the barrier, if there is a component of high excitation 
energy in the n-C motion. Of course the above two mechanisms are not independent but 
mutually related. More analyses are necessary to obtain clear understanding. 

Recently, three-body calculations for the fusion reaction of the same system have been 
reported. These studies, however, report the conclusion opposite to ours, namely the 
strong enhancement of the fusion probability around the barrier energy |6j7]. Let us next 
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consider a reason for this discrepancy. In these three-body calculations, the continuum 
states are discretized in the momentum space. The partial wave expansion of the n-C 
motion is truncated at very low angular momentum, a few units in the %. We believe 
that this truncation is the origin of the discrepancy, causing even an opposite conclusion. 
Figure El shows the fusion probability calculated with l max = 70 (dashed curve) and that 
with l max = 2 (dotted curve). As is seen clearly, the calculation taking only low partial 
waves gives strong enhancement of the fusion probability at low incident energy. However, 
increasing the number of partial waves, this enhancement completely disappears and the 
fusion probability is actually suppressed in the entire energy region. Therefore, it is highly 
probable that the fusion enhancement obtained in Refs. |6|7j is a false account derived from 
the calculations which did not converge with respect to the partial waves. 

We have found that the partial wave expansion up to small /-value gives a qualitatively 
correct answer when one includes only Coulomb breakup mechanism switching-off the 
n-T nuclear potential. Therefore, it is the n-T nuclear potential that does not allow the 
partial wave truncation. A similar behavior is reported in Ref. [12] which points out that 
the inclusion of the high partial waves is necessary to obtain converged breakup cross 
section at low incident energies. 

5. OUTLOOK 

To calculate the fusion cross section, we need to make calculations for reactions at finite 
total angular momenta. This is under progress in the body-fixed frame descrption, where 
the computational cost scales linearly to the total angular momentum. 

For the three-body direct reactions, coupled-channel frameworks discretizing the con- 
tinuum channels in momentum space have been developed and successfully applied. The 
present real-space approach for the three-body problem is more straightforward but com- 
putationally expensive. As the rapid progresses of the computational resources, we expect 
the real-space approach will be more powerful in the future. We would like also to men- 
tion that we have recently developed an alternative real-space approach for the breakup 
reaction employing the absorbing boundary condition|13|. 
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